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Abstract. FIRBACK is a one of the deepest surveys performed at 170 fim with 
ISOPHOT onboard ISO, and is aimed at the study of cosmic far infrared back- 
ground sources. About 300 galax;ies are detected in an area of four square degrees, 
and source counts present a strong slope of 2.2 on an integral "logN-logS" plot, 
which cannot be due to cosmological evolution if no K-correction is present. The 
resolved sources account for less than 10% of the Cosmic Infrared Background 
at 170 ^m. In order to understand the nature of the sources contributing to the 
GIB, and to explain deep source counts at other wavelengths, we have developed 
a phenomenological model, which constrains in a simple way the luminosity func- 
tion evolution with redshift, and fits all the existing deep source counts from the 
mid-infrared to the submillimetre range. 

1 Introduction 



The Cosmic Infrared Background (GIB), due to the accumulation of galaxy 
emission at all redshifts along the line of sight in an instrument beam, is a 
powerful tool for studying galaxy evolution. FIRBACK (|l^, (D), one of the 
deepest surveys performed at 170 fim, is aimed at the study of the GIB, in 
two complementary ways: 

• study the resolved sources (this paper) 

• study the background fluctuations (Q, Q) 

Throughout this paper we use a cosmology with h — 0.65, H = 1 and yl = 0. 
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2 The FIRBACK Survey 

FIRBACK, is a survey of 4 square degrees in 3 high galactic latitude fields, 
chosen to have as low an HI column-density as possible, typically Nh — 
10^°cm~^, and if possible multiwavelength coverage. Observations were car- 
ried with ESA's Infrared Space Observatory (ISO, with the ISOPHOT 
photometer in raster mode (AOT P22) with the C200 camera and C_160 
broadband filter centered at A = 170 /zm. A detailed description of the reduc- 
tion, data processing, and calibration will be discussed in whereas the 
analysis of the complete survey will be discussed in |^ . 

3 Source Counts at 170 fim 

Preliminary FIRBACK integral source counts at 170 (Fig. [^), not cor- 
rected for incompleteness, show a strong slope of 2.2 between 120 and 500 
mjy. This strong slope is not explained by a K-correction or cosmological 
evolution alone: both must be present; the K-correction is the ratio, at a 
given wavelength, of the emitted flux over the redshifted flux. 
Non (or low) evolution scenarii, or extrapolation of IRAS counts, are not 
able to reproduce the observed counts. For illustration, we plot in Fig. |^ the 
non-evolution model from |^ and the evolution model A from ||ll[. On the 
other hand, evolutionary models from |9| and from (model E: evolution 
-I- ULIRGs) give a better agreement. At this wavelength, FIRBACK sources 
account for only 3% of the background. 

4 Modelling the Evolution of Galaxies 
4.1 Method 

Our philosophy is to make a phenomenological model that explains all the 
observed deep source counts and reproduces the CIB in the mid-IR to submil- 
limetre range. One way to do this is to constrain the evolution with redshift 
of the luminosity function (LF) in the infrared, given: 

• templates of galaxy spectra 

• the energy density available at each redshift 

For the first point, we used template galaxy spectra based on IRAS colors 
modified to account for recent ISO observations, in particular the absorption 
feature near 10 /im at high luminosity; PAH features are present in the mid- 
infrared 0, even if their strength seems larger than the observations: this is 
not a problem because the right amount of energy is present in each peak, 
and we convolve the spectrum with the filter spectral response (Fig. bh. For 
the second point, we used the inversion of the CIB spectrum by |10|] and 
[ pO| , which gives with good accuracy the energy density available at redshifts 
between 1 and 3. This energy density is the integral of the luminosity function 
at each redshift, but there is no unique solution for the LF shape. 
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Fig. 1. FIRBACK integral source counts at 170 fim not corrected for incomplete- 
ness. Models from : A (dot) with evolution, E (solid) with evolution + ULIRGs. 
Models from M: without evolution (dot-dash), with evolution (dash). 
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Fig. 2. Template spectra used for simulating source counts; from bottom to top: 
L = 10^, 10^°, 10", 10^^ 10^3 Lq. 
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4.2 Evolving Luminosity Function 

The question is: how to evolve the LF with redshift ? 

Fig. represents the LF of pl| at a redshift z=0, renormahzed toh = 0.65 and 



having the same integral as 1 23 . Our constraint on the LF redshift variation 
is the energy density available at each redshift 0. Usually, authors apply 
a pure density evolution to the LF as a function of redshift (i.e. a vertical 
shift), or a pure luminosity evolution (a horizontal shift). This does not work 
for FIRBACK source counts: There is no alternative other than adding the 
evolution to one part of the LF only. This part is constrained by IR and 
submm observations: this is the bright end of the LF. Fig. ^ represents our 
decomposition of the local LF into two parts: 

• left part: "normal" galaxies 

• right part: ULIRG's, centered on a luminosity LjjLUfa — 2.0 x IO^^Lq, 
where Ljjlirg is the free parameter; we get the same value as 

In our model, "without evolution" means that the local LF (Fig. |a) is taken, 
and is the same at every redshift. In the evolutionary scenario, only the 
ULIRG part moves, in such a way that the integral of the LF equals the 
constraint given by the GIB inversion at each redshift. The maximum is 
reached at z ~ 2.5, and Fig. ||b represents the LF at this redshift. We neglect 
at this stage the evolution of "normal" galaxies, which do not play a crucial 
role in the mid-infrared to submillimetre wavelength range. 




Fig. 3. a: Luminosity Function at z=0 (solid line); normal galaxy (dot-dash); 
ULIRG (dash-dash), b: Luminosity Function at z=2.5 (solid line); normal galaxy 
(dot-dash) ULIRG (dash-dash) and local LF (dots). 
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4.3 Model of Source Counts at 170 fxra 



The model at 170 /Ltm, together with our data, is presented in Fig. ^. The 
brightest point of the observed counts is compatible with all our models, in 
particular the non-evolutionary scenario, which is expected for local sources. 
We also show the effect of the K-correction, which steepens the integral source 
count slope. Our evolutionary scenario fits the data within the error bars. 
Most of the background is expected to be resolved into sources once we are 
able to detect sources at the mJy level. This wavelength region, in which 
evolutionary effects are particularly important and where there are good 
prospects for detecting higher redshift sources because of the K-correction, is 
nowadays probably the best-suited range for probing galaxy evolution from 
space in the far-infrared. 
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Fig. 4. Observed Counts and Models at 170 /im, with (solid line) and without 
(dashed line) evolution. The model with evolution and without K-correction is the 
dot-dash line. 



4.4 Models of Source Counts at other wavelengths 



Data at 15/ito [Q, 90 /im ||] and 850 /zm with our models are presented 
in Fig. |[ Fig. ^ and Fig. |7 respectively. 

At 15 fj,m, both the slopes and the little "waves" in the counts appear in the 
models with the combined effects of the K-correction and the evolution. At 
90 nm, the K-correction does not emphasize the differences between the sce- 
narii of evolution or non-evolution. At 850 fj,m, the "non smooth" appearance 
of our model is due to the discretization of the LF. 
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Fig. 5. Models at 15 /im with (solid line) and without (dashed line) evolution, and 
observed counts from Elbaz et al. |7| (thin lines). 
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Fig. 6. Models at 90 fim with (solid line) and without (dashed line) evolution, and 
observed counts from Efstathiou et al. (squares) and IRAS counts |^] (solid area) 
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Fig. 7. Models at 850 /im with (solid line) and without (dashed line) evolution, and 
observed counts from Barger et al. (diamonds). 

5 Discussion 

5.1 The bright end luminosity function evolution model 

Our model of the evolving LF scenario fits most of the existing deep survey 
data from space (mid and far infrared) and ground (submillimetre) . It is also 
compatible with the observational estimate of the LF in The strong 
observational constrain of multiwavelength source counts is thus explainable 
by a simple evolutionary law of the LF: the Bright End Luminosity Function 
Evolution (BELFE) model. One simple model is in agreement with all up- 
to-date observables and reproduces the background. 

5.2 Redshift and Nature of the sources 

Another crucial observational test for our model is the predicted vs observed 
redshift distribution. Although the statistics are poor, we have some evidence 
that most of the ISOCAM sources lie at redshift between and 1.4 with a 
median at 0.8 and that most of the SCUBA sources lie at redshifts greater 
than 2 Q. Our predicted redshift distributions at these two wavelengths are 
in agreement with the existing observations. 

What about FIRBACK 170 iim sources ? Our predicted redshift distribution 
shows that most of the sources lie at redshifts below 1.5, with a median 
comparable to that of the ISOCAM sources. This means that we are sensitive 
both to local sources and sources beyond redshift z ~ 1. All the FIRBACK 
sources with known redshifts (less than 10) ||], are in this range. Two 
submillimetre sources are at z > 1, and a few visible sources are at z ~ 0.2. 
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It is difficult to address the question diff'erences between PHOT and CAM 
sources, because their redshift distributions (both expected and observed) 
are similar. About half of the 170 /im sources have 15 fim counterparts. 

6 Conclusion 



We have presented results from the FIRBACK survey, one of the largest 
ISO programs, dealing with resolved sources of the GIB: our source counts 
at 170 /im show strong evolution. This evolution is explained by a simple 
law involving the redshift of the luminosity function, which differs from pure 
density or luminosity evolution: the "Bright End Luminosity Function Evo- 
lution" model. The model fits all the existing source counts at 15, 90, 170 
and 850 fim, and also predicts a redshift distribution in agreement with the 
(sometimes sparse) observations. This powerful tool, based on observational 
constraints on the GIB spectrum inversion and the local Luminosity Func- 
tion, and on assumed template galaxy spectra, not only agrees with existing 
data but also is able to make useful predictions on source counts or GIB 
fluctuations. These predictions may be useful for planning the utilization of 
major telescopes of the future, such as SIRTF, Planck, FIRST and ALMA. 
All the FIRBAGK materials are available at: http://wwwfirback.ias.fr. 
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